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I.  INTRODUCTION 


The  complete  history  of  the  experimental  development  of  solid 
mechanics  in  general  and  plasticity  in  particular  is  traced  in 
the  treatise  of  Professor  James  F.  Bell,^*  published  in  1973. 

The  practical  developments  in  plasticity  are  rather  recent,  most 
of  them  occurring  after  1920.  The  first  complete  books  on  plas- 
ticity in  English  were  published  by  Prager  and  Hodge^  and  Hill^ 
as  recently  as  the  early  1950 's  and  to  the  writer's  knowledge, 
the  first  and  only  book  devoted  completely  to  dynamic  plasticity 
was  published  by  Cristescu^  only  eight  years  ago.  The  field  of 
dynamic  plasticity  is  therefore  relatively  new  compared  to  dynam- 
ic elasticity  which  is  several  hundred  years  old. 

The  vulnerability  of  structures  depends  upon  how  much  deformation 
they  can  take  and  still  remain  serviceable.  Vulnerability  implies 
plastic  deformation  and  failure.  Therefore  in  order  to  handle  the 
vulnerability  problem  a knowledge  of  plastic  deformation  of  struc- 
tures is  a prerequisite.  Most  of  the  Army  work  in  vulnerability 
has  of  necessity  been  empirical  because  the  state  of  the  art  in 
plasticity  theory  has  not  kept  up  with  the  necessities  of  prac- 
tical Army  problems.  The  classical  theories  of  plasticity  devel- 
oped over  the  past  twenty  or  thirty  years  up  to  the  early  1970 's 
are  difficult  to  apply  and  leave  something  to  be  desired  in  the 
way  of  accuracy  in  describing  material  behavior.  A summary  of  the 
significant  vulnerability-oriented  calculations  for  simple  struc- 
tures is  contained  in  a recent  paper  by  Westine  and  Baker. ^ It 
is  only  recently  that  Bell^  has  developed  a new  experimentally 
based  general  theory  that  can  be  applied  readily  to  general  large, 
deformation  plasticity  problems  of  the  type  contained  in  vulner- 
ability studies. 

It  is  the  purpose  of  this  report  to  present  the  practical  vulnera- 
bility-oriented approaches  using  plasticity  theory  and  then  intro- 
duce some  calculations  which  are  based  upon  Bell's  new  theory. 


*Superscripts  refer  to  references  listed  at  the  end  of  the  report. 
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II.  ENERGY  THEORIES  IN  VULNERABILITY 


A.  Energy  approaches  in  general 


B. 


The  capacity  of  a structure  to  absorb  energy  can  be  a very  useful 
concept  in  assessing  its  vulnerability.  Three  energy  type  approa- 
ches that  have  led  to  fruitful  results  in  vulnerability  problems 
are  the  variational  approach,  the  asymptotic  approximation,  and 
the  direct  energy  equalization.  All  three  use  the  plastic  energy 
absorbed  in  the  structure  under  a given  configuration  called  the 
failure  configuration  or  the  mode  of  failure.  This  energy  absorbed 
will  be  denoted  by  V,  where  V is  a function  of  the  geometry  of  the 
structure,  its  material  properties  and  the  deflection  component 
in  the  direction, of  the  load.  These  three  approaches  will  be  dis- 
cussed in  the  order  of  increasing  complication.  The  objective  of 
all  the  approaches  is  to  compute  the  final  plastic  deflection  dis- 
tribution of  the  structure  as  a function  of  an  impulse,  load  or 
energy  magnitude.  Knowing  this  deflection  we  can  then  be  in  a 
position  to  judge  whether  or  not  this  deformation  will  constitute 
a failure. 


Direct  energy  equalization 

In  this  approach  the  energy  absorbed  by  the  structure  is  equated 
to  the  energy  directed  toward  the  structure  either  from  an  explo- 
sion or  from  a series  of  fragments.  The  main  difficulty  with 
this  approach  (once  V is  known  for  a given  failure  configuration) 
is  the  estimation  of  the  energy  given  to  the  target  by  the  explo- 
sion or  fragments.  For  blast,  in  several  past  reports®'^  the 
writer  has  taken  this  explosive  energy  to  be  the  energy  flux 
(energy  per  unit  area)  in  the  explosive  wave  multiplied  by  the  pro- 
jected area  of  the  target  facing  the  explosion.  For  fragments, 
assuming  that  perforation  takes  place,  the  energy  absorbed  by  the 
structure  is  computed  by  subtracting  the  residual  kinetic  energy 
of  the  broken  up  fragments  after  perforation  from  the  initial  kinet 
ic  energy  of  the  approaching  fragments. 


Since  this  estimation  of  blast  energy  does  not  include  any  inter- 
action effects  between  target  and  blast  wave  it  is  therefore  bound 
to  be  in  error.  Nevertheless  it  has  led  to  a basic  understanding 
of  the  role  of  energy  in  predicting  the  form  of  the  isodamage  cur- 
ves. At  the  present  state  of  the  art,  we  are  at  the  stage  where 
this  approach  can  and  should  be  considerably  refined. 


C.  Asymptotic  approximation 

The  asymptotic  energy  approach  has  been  formulated  by  Westine  and 
Baker^  in  a very  recent  report  and  can  be  considered  one  step  more 
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complicated  than  the  direct  energy  equalization  approach  discussed 
above.  They  equate  the  energy  absorbed  in  the  structure  to  two 
separate  characteristic  energies  of  the  loading  in  order  to  narrow 
down  on  the  damage  characteristics.  The  first  of  these  character- 
istic energies  is  the  kinetic  energy  imparted  to  the  structure  for 
short  duration  impulsive  loads.  The  second  characteristic  energy 
is  the  work  performed  by  the  peak  load  moving  through  the  distance 
that  the  structure  deforms.  Thus  two  extremes  or  asymptotes  for 
the  energy  imparted  to  the  structure  are  calculated.  One  of  these 
is  in  the  very  short  time  impulsive  loading  regime  and  one  is  in 
the  long  duration  quasi  static  regime.  In  this  way  the  damage 
characteristics  are  approached  from  both  ends  of  the  loading  spec- 
trum. 

D.  The  variational  approach 

In  the  energy  equalization  approach  no  account  was  taken  of  the 
interaction  of  the  blast  load  with  the  structure  so  loading  dis- 
tribution and  timewise  effects  were  neglected,  in  the  asymptotic 
approach  some  account  was  taken  of  the  load  distribution  but  time- 
wise  effects  were  only  considered  in  a limiting  way.  In  the  var- 
iational approach  both  timewise  effects  and  load  distribution 
characteristics  are  included.  The  variational  equations  for  a 
structure  with  a given  failure  configuration  were  derived  in  a 
previous  report. For  a given  failure  configuration  the  system 
reduces  to  a single  equation  for  the  lateral  deflection  vhich  is 
composed  of  an  inertia  term  depending  upon  the  mass  distribution, 
a stiffness  term  depending  upon  the  plastic  plus  elastic  energy 
absorbed  in  the  given  configuration  and  a loading  term  which  is 
dependent  upon  the  pressure  distribution  over  the  structure.  The 
equation,  in  general,  is  nonlinear  in  the  deflection,  therefore 
this  approach,  even  though  more  accurate  in  principle  than  the 
ones  described  in  Sections  B and  C,  can  be  considerably  more  com- 
plicated if  higher  order  terms  in  the  deflection  are  involved. 

There  are  a number  of  practical  cases,  however,  which  only  involve 
up  to  quadratic  terms  in  the  deflection  and  these  will  be  discussed 
later  in  this  report. 


III.  MATHEMATICAL  DETAILS  OF  THE  ENERGY  APPROACHES 


A.  Direct  energy  equalization 


1.  General  equations 

7 

The  writer  has  proposed  using  the  following  energy  equation 
for  predicting  the  isodamage  characteristics  of  structures: 


where 


P I 


Ef  = 


+ 


3^2 


n V 


2 M 
A 


[11 


energy  density  (energy  per  unit  area) 
which  is  absorbed  by  the  target 
structure  in  undergoing  a certain 
degree  of  damage 
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P = side  on  pressure  in  the  blast  impinging  on 
the  structure 

1 = side  on  impulse  in  the  blast  impinging  on  the 

structure 

p = mass  density  of  the  ambient  air 
o 

c = sound  velocity  in  the  ambient  air 
o 

n = number  of  fragments  hitting  the  target 
V = velocity  of  the  fragments  at  the  target 
M = mass  of  each  fragment  which  hits  the  target 
A = effective  area  over  which  the  damage  takes  place 
If  we  multiply  equation  [1]  by  A we  obtain 

^total  ~2pc^  t^nVM 

o o 

^total  \last  ^fragments 


The  above  relations  state  that  the  total  energy  absorbed  by  the  struc- 
ture undergoing  a given  level  of  damage  is  equal  to  tne  blast  energy 
plus  the  energy  imparted  by  the  fragments.  The  rationale  and  physical 
principles  behind  this  relation  are  discussed  thoroughly  in  the  ear- 
lier reference.^  If  there  are  vaporific  effects,  i.e.  if  the  frag- 
ments perforate  the  structure,  hit  an  internal  component  or  another 
piece  of  structure  and  produce  internal  pressure  because  of  melting 
and  vaporization,  equation  [2]  still  holds  because  it  is  tacitly 
assumed  that  all  energy  from  the  blast  or  fragments  is  absorbed  by 
the  structure  and  vaporific  effects  just  represent  a conversion  of 
fragment  kinetic  energy  into  heat. 

2.  The  experimental  approach 

Relation  [2]  can  be  applied  in  two  ways.  The  first  way  is  to  use  an 
experimental  approach.  A test  must  be  performed  on  the  structure  to 
obtain  one  point  on  the  damage  curve,  i.e-  perform  either  a blast  or 
fragment  test  to  obtain  a given  level  of  damage  on  a structure.  If  a 
pure  fragment  test  is  performed,  then  the  result  can  be  plotted  as  a 
point  on  the  V-M  curve  as  shown  in  Figure  1. 
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Pig.  1 Fragment  isodamage  Curve  with  No  Blast  Effects 

With  n fragments  each  of  mass  M hitting  the  target  at. velocity  V the 
energy  imparted  to  the  target  will  be 

E = h nMV^  [4] 


So  the  isodamage  curve  is  described  by  the  equation 

This  isodamage  curve  is  shown  as  the  curve  in  Figure  1 drawn  through 
the  experimental  point.  The  value  of  E given  by  equation [4]  is  taken 
as  the  total  energy  necessary  to  do  damage  of  the  given  level.  If 
blast  effects  were  present  then  we  could  describe  the  blast  energy 
as  • Then  the  fragment  energy  necessary  to  afflict 

the  "same  damage  level"*  would  be 


E = E - 

^f  ® 2p  oCo 


A 


[6] 


The  isodamage  curve  for  the  fragments  would  then  be  described  by 


V = 


Ef 


[7] 


h n M. 


* This  has  been  put  in  quotes  since  fragment  damage  and  blast  damage 
do  not  always  afflict  the  same  type  of  damage. 
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similarly  if  a blast  test  were  run  on  the  structure  to  do  damage  at 
a given  level,  this  point  could  be  plotted  on  a blast  isodamage  curve 
as  shown  in  Figure  2 below. 


Fig.  2 Blast  Isodamage  Curve  with  No  Fragment  Effects 
The  energy  for  doing  damage  at  the  given  level  is  given  by 

P I A 


E = 


^ P 

o 


[8] 


So  the  isodamage  curve  is  given  by  the  relation 

E 


[9] 


P = 


2 P-  c. 


This  curve  is  shown  as  the  solid  line  in  Figure  2.  Now  if  fragments 
were  present  their  energy  would  be  described  by  ^nMV^  and  the  blast 
energy  necessary  to  do  the  same  damage  would  be 

,2 


Eg  = E - ^5  n MV^ 


[10] 


The  new  isodamage  curve  for  blast  would  then  be  given  by  the  relation 
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p 


I 


[11] 


B 


3.  The  theoretical  approach 


A 


The  second  way  to  apply  the  basic  equations  is  to  calculate  the 
total  energy  under  a given  failure  pattern  at  a given  damage  level. 
Several  earlier  references®'^  contain  methods  for  computing  blast 
damage,  the  most  complete  presentation  being  given  in  Reference  9. 
Fragment  damage  and  blast  damage  occurring  after  fragment  damage 
is  discussed  in  Reference  11.  The  energy  necessary  to  do  frag- 
ment damage  can  be  computed  from  the  equation 

E = H n M 

xn 

Eq.  [121  is  a perforation  equation  in  which  n is  the  num- 
ber of  fragments  hitting  the  target,  M is  the  fragment  mass  and 
Vxn  is  the  ballistic  limit  velocity  (i.e.  the  minimum  velocity 
necessary  to  perforate  the  target  - see  Ref.  12) . The  energy  nec- 
essary to  do  a certain  blast  damage  is  given  in  Reference  9 for 
various  patterns  of  failure.  The  pattern  of  failure  must  be  assum- 
ed in  advance. 


4.  Checking  the  theory 

The  formulas  in  the  previous  sections  can  be  validated  by  using 
existing  data  on  damage  of  structures.  The  way  to  achieve  this 
is  to  construct  the  isodamage  curve  by  using  a single  data  point 
from  an  existing  test  and  then  checking  to  see  if  other  points 
fall  along  the  same  curve.  To  check  the  theoretical  approach  of 
computing  the  energy  absorbed,  the  energy  can  be  computed  by  assum- 
ing a given  pattern  of  failure,  constructing  the  isodamage  curve 
from  this,  and  then  checking  experimentally  to  see  if  the  pressure 
and  impulse  or  velocity  and  mass  values  fall  along  this  curve  for 
equivalent  damage. 

B.  Asymptotic  approximation 

As  an  alternative  to  solving  the  details  of  each  problem,  Westine  and 
Baker^  have  developed  a procedure  using  the  absorbed  plastic  energy 
to  compute  the  asymptotes  for  the  impulsive  and  quasi-static  regions. 

They  determine  these  asymptotes  by  equating  the  internal  plastic  work 
(or  strain  energy)  first  to  the  kinetic  energy  imparted  to  the  struc- 
ture to  get  the  asymptote  for  the  impulsive  loading  regime  and  then 
to  the  work  performed  by  the  peak  force  to  get  the  asymptote  for  the 
quasi-static  regime.  Let  V be  the  plastic  energy  absorbed  (or  inter- 
nal work) , This  value  of  V is  given  for  cylindrical  shells  by  eq. [82] 
and  Figs,  22-27  of  the  writer *s  1970  report  and^£or  lifting  surfaces 
®9^-[44]  - [46]  and  Table  2 of  his  1971  report  . Equating  this  energy 
absorbed, V, to  the  kinetic  energy  imparted  to  the  structure  for  the 
impulsive  loading  .regime,  the^ impulse  per  unit  area,  I,  is  given  in  terms 
of  the  energy  by  the  relation 
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I 


/ 


.UL— 

f (A)dA 


[13] 


where  f(A)  is  the  impulse  distribution  on  the  structure  and  jj. 

is  the  mass  per  unit  area  of  the  structure. 


[14] 


[15] 


In  the  quasi-static  loading  regime  the  work  done  by  the  peak  load  is 

W = j'  P (A)  w(A)  dA 
A 

Equating  this  to  V we  obtain  for  the  quasi-static  region 

V = W 

(V  is  A function  of  the  deflection  w) 

where  P(A)  is  the  spatial  load  distribution  and  w is  the  lateral 
deflection.  Some  special  cases  for  both  axi symmetric  and  nonaxisym- 
metric  collapse  of  cylindrical  shells  were  considered  by  the  writer 
some  years  ago. 

C.  Variational  approach  to  the  blast  problem 
1.  General  equations 

There  is  another  approach  to  the  damage  problem  which  looks  closer 
at  the  individual  structure  and  follows  the  damage  mechanism  as  it 
occurs.  This  approach  can  best  be  illustrated  by  an  example.  Con- 
sider a shell  subjected  to  an  enveloping  blast.  By  using  varia- 
tional principles  as  given  in  an  earlier  report^^  the  equation  for 
the  plastic  radial  deflection,  w,  cylindrical  shell  under  non- 

axisymmetric  loading  can  be  written 

[16] 


w 


= w_(t)  f (A) 

w 


where  f (A)  is  the  spatial  distribution  of  def lectioj^over  surface  A and 
w^  (t) , ¥he  timewise  part  of  the  deflection,  satifies  [17] 


o 

w 


r 


j 


jjL  (A)  dA 


w 


+ 


8v 

8 w 


P(A,t)  f (A)  dA 
w 


A A 

or  written  in  more  familiar  single  degree  of.  freedom  notation 

m X t R (x)  = P (t) 

e e ^ , 

where  has  been  replaced  by  x and 


[18] 


m 


= j ^ 


T (A)  dA  = the  generalized  mass 
w 

= distribution  of  deflection  w over  sur- 
- face  A 

(x)  = = the  resistance  (elastic  or  plastic)-  it 


19] 


can  be  approximated  by  a power  series  in  x 
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and 


where 


P„(t)  = 

e 


P(A,t)  f (A)  dA 
w 


[20] 


P(A,t)  = Po(A)f(t) 

Pq  (A)  = spatial  load  distribution 
f(t.)  = timewise  load  distribution 

IJi  = mass  per  unit  area  of  structure 
2,  Form  of  the  solution  for  a long  shell  (see  Ref.  10) 

For  a perfectly  plastic  material  in  which  1 


[21] 


= deflection) 


L = shell  length ,^Wj. 

O A 

where  h = shell  thickness,  CJ5  = yield  stress  in  pure  tension. 


(a  = shell  radius, 

3y 

0w 


f (A)  dA 
w 


[22] 


If 


we  further  limit  the  discussion  to  exponential  timewise  loading 
(i . e.  P (A,  t)  =p  (A)  f (t)  =p  f (A)  e~^' ),  then  the  equation  of  mo- 
tion [18]  takes  the  same  ^form  as  eq 

t/T 

P e” 


and  Baker , ^ i . e . 
where 


[1]  of  the  paper  by  Westine 
f = m X [23] 


P = 


P^(A) 

o 


f (A)  dA 
w 


o p 


f (A)  dA 
w 


[24] 


= a h 
s 


^7^ 


f (A)  dA 
w 


^ fp (A) =spatial  load  distribution 

' j M 

A 

5 


m 


The  solution  curve  is  exactly  the  one  shown  by  Westine  and  Baker  as 
Figure  2 of  their  paper  with  the  above  parameters  which  depend  upon 
the  deformed  shape,  spatial  load  distribution,  fp(A) 

and  the  other  physical  cnaracteristics  of  the  shell.  This  curve  is 
shown  in  Fig.  3 
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In  this  curve 


I 

X 

3.  Interpretation  in  terms  of  conventional  P-I  isodamage  curves 

An  interpretation  of  this  scaled  P-I  curve  as  it  relates  to  damage 
problems  is  certainly  in  order  here.  This  curve  (fig.  3)  really 
gives  the  solution  to  the  plastic  problem  in  nondimensional  form. 
For  a given  set  of  values  of  P,  f,  I,  m we  can  obtain  the  value  of 
I^/Xmf  resulting  from  a certain  value  of  P/f  and  then  calculate  the 
value  of  the  maximum  deflection,  X,  from  this  result.  The  conven- 
tional isodamage  curve  is  a plot  of  pressure  vs.  impulse  as  shown 
in  Figure  4. 


OO 

; 


op  . 

P e-  dt  = J P„  e"  at  J f (A)  £„(A)  dA 


0 0 
= maximum  deflection  (i.e.  (wQ)max.) 


Pressure  ; 

i 

P ' 

O i 


Damage  greater  than  X 


Damage  less  than  X 


OO 


Impulse  I 


. r p 

J O 


- t/T 


dt  = 


P T 


Fig.  4 Conventional  Isodamage  Curve 

The  conventional  curve  of  Figure  4 is  obtained  from  the  nondimen- 
sional solution  curve  of  Figure  3 as  follows: 

a.  Choose  values  of  P/f  from  Fig.  3 and  calculate  Pq  from  the 
relations 
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= f (p/f) 


and 


f f (A)  f^(A)  dA  [26] 


b.  Pick  off  values  of  I^/Xmf  from  Fig.  3 corresponding  to  the  partic- 
ular values  of  P/f  selected  in  "a'  above.  Compute  l from  the  relation 


X m f 


Value  obtained  from  curve  = l [27] 


Assume  X for  a given  damage  level,  thus  o _ 

I = I X m f 

[ f f (A)  f (A)  dA  ]^  = i X m f 

o J p w 


therefore 


I 

o 


^ i X m f 

r f (A)  f (A)  dA 
p w'  ' 


[28] 


Note  that  the  conventional  isodamage  curve  calculated  in  this  way 
(as  contrasted  with  that  shown  in  Section  IIIA)  will  not  have  zero 
asymptotes,  i.e.  it  will  give  infinite  pressure  for  some  finite  value 
of  impulse  instead  of  infinite  pressure  for  zero  impulse^  Likewise  it 
will  give  infinite  impulse  for  some  finite  value  of  pressure  instead 
of  infinite  impulse  for  zero  pressure.  This  is  therefore  a more  accu- 
rate way  to  obtain  the  isodamage  curve  for  blast  on  a given  structure 
if  the  load  distribution  is  known. 


IV.  VARIATIONAL  SOLUTIONS  FOR  BEAMS  AND  PLATES 


A.  General  background 


The  general  form  of  the  variational  equation  is  given  by  equation 
[18].  If  the  resistance  function,  is  a constant  then  the  equa~ 
tion  takes  the  same  form  as  the  simple  Westine-Baker^  equation 
(i.e.  eq.  [23])-  For  these  cases  the  curve  of  Westine  and  Baker 
(i.e.  Fig.  3)  can  be  used  as  a nondimens ional  isodamage  curve  as 
long  as  the  proper  interpretation  is  given  to  each  of  the  terms. 


B.  Plastic  deformation  of  cantilever  beams 

c 

Westine  and  Baker  have  computed  the  plastic  strain  energy  of  a 
cantilever  beam.  This  energy,  V,  is 


7T  O'  b h^  w 
..,.X o 

16  L 


[29] 
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where 


Oy=  yield  stress 

b = width  of  beam 

h = thickness  of  beam 

L = length  of  beam 

w = maximum  tip  deflection 
o 


Using  a deflection  shape  of  w = w ( i - qq3  7T x 


2 L 


The  resistance  function,  for  this  case  turns  out  to  be 

e Ow  16  L 

The  generalized  mass,  m ° is 


[30] 

[31] 


m 


/X  f^(A)  ciA 


e 

= b 


Ip  h ( 1 - cos  ^ ^ dx  [32] 

r 2 L 


^ Q = .23  b h L 

For  a uniformly  distributed  load  the  generalized  loading  function 
is 


P = 


f (A)  f (A)  dA 


w 


fp(A)  = 1 


[331 


Thus 
P/f  = 

and 


P (.36  b L ) 16  L 
° -2 


TTCTy  b h 


X m f 


a 

Y 


2 2 
(.36  b L ) 


P = P (.36  b L 
o 


16  L (.36) 
-2 


7Th‘ 


a 

Y 


X (.23  bp  h L)  rra^  b h' 


p 


h w 


16  L 


[34; 

(1.83) 


[35] 


L 


(2.88) 


So  the  curve  of  P/f  vs.  1‘^/Xmf  converts  for  a cantilever  to  a curve 


of 


1.83 


vs 


2.88 


Y 


ph  (J 


cantilevers  made  from  6061-T6  aluminum  which  were  12  inches  long 
and  .051  inches  thick  were  tested  at  BRL.^^  For  this  material 
CTy  = 40,000  psi,  P = .000255  lb.  sec.Vin.'^.  The 

asymptotes  of  the  isodamage  curve  are  calculated  directly  as 
follows: 
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For  the  impulsive  loading  regime 


/ X m f 


2.0 


so 


2.88 


2.0 


[36] 


h w O' 
o y 


therefore 


w 


1.44 


L 

h 


IP 


y 


For  the  quasi-static  loading  regime 


so 


1.83 


L, 

T 


= 1 


or 


P / f = 1 
o 


. 54 


[37] 


The  value  for  the  quasi-static  regime  checks  with  the  Westine- 
Baker  value^  since  they  take  account  of  the  deflection  distribution 
in  their  calculation  of  the  quasi-static  asymptote.  The  value  cal- 
culated for  the  impulsive  regime  does  not  check  since  they  omitted 
the  effect  of  deflection  distribution  on  the  Kinetic  energy.  If 
the  energy  principle  is  derived  from  the  variational  equation,  it 
is  found  that^^  x 

38] 

p (t)  dX  - . T.  ^ J 


m 

R (x)  dx 
e 


where 


P (t) 
e 


■ / 


0 


P (A,  t)  f (A)  dA 
w 


R^(x)  = 


dv 


[39] 


Under  thej^  circumstances  the  result  for  the  impulsive  loading 


region  is 


where 


H 


me 


2 m 


[40] 


= strain  energy  absorbed  by  the  structure  (inter- 
nal work  done  by  the  structure) 

T 


H 


me 


m 


/ 

/ 


Pg(t)  dt  ( generalized  impulse) 

2 ^41] 

f (A)  dA  , ^ . 

' ^ V the  generalized  mass/ 
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The  nondimens ional  isodamage  curve  for  the  cantilever  is  then  found 
by  replacing  the  ordinate  and  abscissa  of  Figure  3 by 


1.83 


2 

and  2 . 88  ^^o 

P *’  "o  O-y 


respectively 


The  conventional  isodamage  curve  for  this  cantilever  is  shown  in 
Fig.  6. 
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There  are  several  basic  principles  which  can  be  learned  from  this 
example.  Firstly,  the  isodamage  asymptote  for  pressure  is  a con- 
stant for  the  perfectly  plastic  beam,  i.e.  it  is  independent  of 
the  level  of  damage.  Secondly,  the  asymptote  for  impulse  is  de- 
pendent upon  the  damage.  The  pressure  asymptote  really  defines  a 
collapse  pressure  for  the  beam  under  static  loading. 

C.  Plastic  deformation  of  uniform  simply  supported  beams 


For  the  case  of  a perfectly  plastic  simply  supported  beam  of  rec- 
tangular section  the  strain  energy  is  given  by^ 

8 M w 


V 


where 


M 


O'  b h L 

L yield  moment 


o 


b = width  of  beam 
h = thickness  of  beam 
L = length  of  beam 


Thus 


f = R 

e 


3v 


8 M 


w 


[431 


[44] 


[45] 


TTx 


The  generalised  mass  is  (using  a deflection  shape  of  ^ 

) r 2 r ^ • 2 TTx  - 

I p h sin  — dx 
J [4bj 


w = w sin 

o L 


m 


d V 


^ f^(A)  dA  = b 


.5 


P 


h b L 


For  a uniformly  distributed  load  the  generalized  force  is 


P = P 


f (A)  f (A)  dA 
P ^ 


fp(M  - 1 


[471 


So  P = .64  b L P 


So  for  this  case 


P/f  = ‘’o 


P L 
o 


8 cr  b h 

I 


(.32) 


(J  ^2 
y h 


^2  . ^ (.64  b L)^ 

I / X m f = o 


[481 
^ 2 


w (.SbLnh)  8 (J  bh 

O P ^ ) 


(.64) 


nh  w.  (J  -2 
r o y h 


4 L 


So  the  curve  of 


P/f  vs  I / X m f 

2 


supported  beam  to  a curve  of  p L 

.32  o 


y h 


converts'  for  a simply 
41 

VS  .41  o 


nh  w (J  ,2 
r o V h 
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The  asymptote  for  the  impulsive  loading  region  is 


.41 


So 


w 


n h w (J 
r o ^ y 


.2  { 


) 


L 

IT 


[49] 


h-  “ ’’•V^ 

In  the  paper  by  Westine  and  Baker^  their  2 £ = L that  the 

above  value  comes  very  close  to  their  value.  The  slight  difference 
is  probably  due  to  the  fact  that  the  shape  was  considered  here  in 
computing  the  kinetic  energy.  The  asymptote  for  the  quasi-static 
regime  is  given  by 


.32 


■or 


3 (h/L) 


[501 


which  is  exactly  the  value  given  by  westine  and  Baker. 


D.  Plastic  deformation  of  uniform  simply  supported  plates 

For  the  case  of  simply  supported  plates'  Westine  and  Baker^  find 
that  the  strain  energy  is  given  by  (assuming  a perfectly  plastic 
material  and  plate  dimensions  of  a,  b,  h - width,  length, thickness 


respect Ivelyl . 

V = CJ  h^  w ( — 
y o a 


- ) 
b ' 

7T‘ 


4 


cr  h^  w 

y o 


8 


0“  h w 

y o 


[51] 

5 ) + A cr  h w^ 

h ^3  y o 


The  terms  linear  in  represent  bending  terms  and  those  quadratic 
in  arise  from  tension  in  the  middle  surface  of  the  plate. 

Let  w be  given  by  _ . 7T  x 7T  y 

b [52] 


w = 


w 


sin 


sin 


The  terms  in  the  variational  equation  are  then 


m 


= m 


P = 


u 

/ / 

0 0 

= P 


P 


7T 


TT 


Sin 


a c 

/ f 


X 


sin 


sin 


sin 


dx  dy 


_ p h a b 
= 4 
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dx  dy  = 


a b 


m w 
e o 


0 0 

The  variational  equation  becomes 

4-  ~ so 


O 


m w 


f , w 
1 o 


Again  using  the  typical  exponentially  decaying  pulse 

- , ^ -t/T 

f(t)  = e 


We  obtain 


m w 


+ 


f w 
1 o 


= P e 


-t/T 


- f 


(54] 

P f(t) 

[55] 
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where 


b a 

a b 


O',  h ( - + ^ ) + (T  h 


cr  ( - + ^ ) 

y a b 


4 ^1.2 

+ ^ h 

^3  y 


[571 


This  is  a linear  differential  equation,  the  solution  of  which  is 
as  follows: 


The  initial  conditions  are 
Pi  = V^l/m 


( 0 ) = w^ ( 0 ) = 0 


Thus 


'^o  " m p 


7“ 

P e sin  p,  {t  -T)  dT 


[581 


[59] 


Integrating,  we  obtain 


_ — sin  p,  (t  - T)dT 

m p I 


P T 

w (t)  = 

o'  m p, 


sin  p^t  - p^T  cos  Pj^t  + p^T  e 


-t/T 


[60] 


1 + (p^T) 


'1  ( 1 - cos  p^t  ) 


m p^ 

The  maximum  deflection,  X is  determined  from  the  criterion  that 


i . e . when 


w =0  when  w = (w  ) ( = X) 

o o o max 


P „ , cos  p.t  + p,T  sin  p,t  - e 

£ Pj^T  [ 1 11 1 


-t/T 


[61] 


1 + (p^T) 


sin  Pj^t  = 0 


The  lowest  value  of  t/T  which  satisfies  this  equation  gives  the 
time  at  which  the  maximum  deflection  will  occur.  If  we  let  a=Pj^T 
then  the  above  equation  becomes  -t/T 

cos  a t/T  + CL  sin  gt/T e ^ 

f 1 ^ .2  [62] 


sin  a t/T  = 0 


where 


p^T 


63] 


The  minimum  values  of  t/T  which  satisfy  this  equation  are  plotted 
as  a function  of  p/f  for  various-  CL  in  Figs.  7,8.  The  deflection 
can  then  be  written  as 


X m g 
2 

P T 


[ 


sin  g t/T  - g cos  g t/T 


-t/T 

_a_e ^ [64] 
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P g 

f 


( 1 - cos  a t/T  ) 


In  order  to  compare  the  results  with  the  Westine-Baker  results 
(for  Q.  = 0)  we  plot  P/f  as  ordinate  and  i^/Xmf 

as  abscissa.  The  resulting  curves  are  shown  in  Figure  9a  for 
various  values  of  d , Note  that  for  CT  ^ 0 ^ 

there  is  no  vertical  asymptote  as  found  by  westine  and  Baker  for 
CC  = 0 ' horizontal  asymptote  for  nondimens ional 

pressure  remains  the  same  for  all  CC  . Figure  9b  contains  the 
same  solution  curves  plotted  with  a different  abcissa . This  figure 
illustrates  how  increasing  CC  increases  the  stiffness  of  the  plate. 
The  parameter  d is  a measure  of  the  tension  in  the  middle  surface 
of  the  plate.  The  solution  curves  are  general  and  can  be  applied  to 
any  system  in  which  the  energy  can  be  expressed  as  a quadratic  function 
of  the  deflection. 


V.  APPLICATIONS  USING  THE  BELL  STRESS- STRAIN  LAWS 


A,  General  background 


For  about  a decade  Professor  James  F.  Bell  of  the  Johns  Hopkins 
University  has  been  developing  techniques  and  producing  extensive 
experimental  stress-strain  data  on  many  materials  under  uniaxial 
loading  for  both  the  static  case  and  the  dynamic  case  involving 
high  rates  of  strain.  Recently  he  has  also  made  extensive  studies 
on  biaxial  loading  in  the  static  regime.  Most  of  Bell's  work  ^ ^ 
which  is  applicable  here  is  contained  in  two  recent  references.  ' 
He  has  developed  both  deformation  and  flow  laws  for  alloys 
under  uniaxial, shear , and  biaxial  loading.  in  the  present  report 
the  writer  will  illustrate  the  application  of  Bell’s  Laws  to  both 
one  and  two  dimensional  cases  of  structures  under  impulse  loading. 


The  one  dimensional  tension-compression  stress- strain  law  given  by 
Bell^  for  plastic  deformation  is 

(J  - E € if  € < €y 

I [651 

C7  = CTy  -h  C v€  €y  ^ 

O'  = oje  - ^ if  € > €c 

An 


In  the  above  equations 

CT  = tensile  or  compressive  stress 

€ = tensile  or  compressive  strain 

^ ”=  elastic  modulus 

= yield  stress 
€y  = yield  strain 

= critical  strain  (given  by  Bell^) 
C - 0 given  by  Bell^) 

parameters  given  by  Bell^ 
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Fig.  9 Nondimensional  Solution  Curves  for  Plates  under 
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Each  of  the  parameters  given  by  Bell  is  associated  with  a partic- 
ular physical  phenomena  of  the  material  as  explained  in  Bell's 
previous  work.  There  is  an  analogous  one  dimensional  shear  law 
which  follows  the  same  form  as  the  one  dimensional  tension-compres- 
sion relation. 

For  biaxial  stress  Bell^  gives  as  the  deformation  law 

T = E r If  r < Ty 

, [661 

T = Ty  f syr-Ty  if  TyC  r < Tc 

T = c Vr  - n f & if  r > Tc 

An 

where  5 = o 

where 
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= octahedral  shear  stress 


= octahedral  shear  strain 

All  ghe  parameters  involved  in  the  above  equations  are  given  by 
Bell  and  explained  thoroughly  in  his  treatise^  and  his  many  earlier 
papers. 

B.  Application  to  one  dimensional  problems-bending  of  beams 
1.  Moment  curvature  relationship 


T 

r 


. JiV«i . 4 . 


X + OV  - Ox  CTy  + 3'^Y 


€x  €y  + 


7x% 

'Z~ 


In  the  classical  theory  of  elasticity  the  bending  moment,  M, 
in  a beam  is  given  in  terms  of  the  curvature,  k of  the  beam 
center  line  by  the  relation 

M = E I K 


where  E is  the  modulus  of  elasticity  of  the  beam  material  and 
I is  the  area  moment  of  inertia  around  the  neutral  axis.  This 
relation  is  derived  from  the  basic  assumption  that  plane  sec- 
tions remain  plane  after  deformation,  i.e.  that  the  strain  is 
given  by^^  ^2 


‘X  - ^ K = Z 

where  Z is  the  distance  from  the  center  line  to  any  fiber 


[691 


within  the  depth  of  the  beam . 
given  by  j 

M = / 


Thus  the  bending  moment  is 
b(x,z)  cr^  Z dz 


'A 

where  b(x,z)  is  the  width  of  the  beam  (which  could  be  a func- 
tion of  the  depth  (z)  and  the  longitudinal  dimension  (x) . For 
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a rectangular  beam  of  constant  section  (width  2B  and  depth  2H) 


r+H 

but 

M = 

2 B 

/ (7  z dz 

Thus 
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X 

3 02w 

E 

€ = E z “ 

X 

M = 
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E B H ^ 2^ 

O X 

E I ^ 

a X 
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where  l = j B 
16 

In  the  theory  of  perfectly  plastic  solids  the  perfectly  plastic 
moment,  is  defined  as  the  bending  moment  obtained  when  the 

cross  section  of  the  is  entirely  in  the  plastic  region,  thus 

for  a rectangular  beam"^  H 

[741 


M ~ 4 B 

o 


[ 


z dz  = 2 B (J  H 

Y 


0 


where  Co  is  the  yield  stress  of  the  material  in  pure  tension. 
The  moment  curvature  relation  for  an  elastic-perfectly  plastic  solid 
rectangular  beam  will  then  be  as  shown  in  Figure  10. 

If  the  Bell  theory  is  employed  then  we  must  concern  ourselves  with 
the  three  regions  enumerated  by  Bell.^  For  a rectangular  beam  it 
is  found  that  ^2 


by  letting  K - 
For  € < €y  H 
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M = 4 B 


z E K dz 


So 


and  €=  KH,  the  strain  at  the 
outer  fiber  z=H 


_M _ 4 E £ 


B H 


b.  For  €y,^  € €c 
M = 4 B 


f- 


!y 

K 

Cy 

K 


dz 


H 


f 


K 


z [ CT  + C 


[751 


( Kz  - €^)  '] 


So 


M 


B H cr. 


4 [ 


J + j ( 1 - ) 


2C(2£  +36  )(€-€„) 

+ y iL_ 


3/2 


15  € (J. 


-31- 


c.  For  € ) €c 
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The  moment  curvature  relation  as  predicted  by  the  elastic-perfectly 
plastic  theory  is  compared  to  that  predicted  by  the  Bell  theory  for 
■6061-T6  aluminum  in  Figure  10.  Note  the  smooth  transition  predicted 
by  the  Bell  theory  from  the  elastic  through  the  elastic-plastic  re- 
gions and  finally  into  the  completely  plastic  regime. 

Calculation  of  energy  absorbed  in  elastic-plastic  deformation  of  beams 


The  elastic  energy  absorbed  in  the  bending  deformation  of  beams  can 
be  written  as^  ^ 

(761 

where  M is  the  elastic  moment,  E the  elastic  modulus,  I the  moment 
of  inertia  and  L the  beam  length.  For  a rigid  plastic  material  the 
plastic  energy  absorbed  can  be  written^ 

u = M G>"w„  dx  r77, 

P Jo  ° 

For  an  elastic-plastic  material  such  as  dealt  with  in  the  Bell  Law 

■ 17 

the  energy  absorbed  xs  given  by 


u 


but 


so 


= L'  f 


►m 


O'  d e ] dv 
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(zK) 


z K 
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r r iTiax 

j dx  / dK  4b  / 
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CT  z dz 


2 dK 
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L (zK) 

= J"  dx  M 


[781 

[791 

max 

M(K)  dK  [801 


0 


0 


0 


0 


If  we  consider  the  case  of  a simply  supported  beam  (i.e.  one  where 
the  bending  moment  and  deflection  are  both  zero  at  the  ends  of  the 
beam)  and  assume  a deflection  pattern  as  follows; 

16  X 

^ [811 


w = 


w 


3 2 3 

( L - 2 L X + X ) 


( this  corresponds  to  the  deflection  pattern  under  a uniform  load) 
then  the  energy  can  be  computed  as  a function  of  Wq/L  (Wq  being  the 
maximum  center  deflection  of  the  beam  and  L being  the  length  of  the 
beam) . 

Figure  11  gives  the  results  of  this  energy  calculation,  A sixteen 
point  Gauss  Quadrature  was  used  to  calculate  both  integrals  shown  in 
eq.  [80] . The  computer  program  that  was  used  to  make  these  computa- 
tions is  shown  in  Appendix  I of  this  report.  Note  that  for  small 
elastic  deformations  U varies  as  (w^/L)  whereas  for  the  larger 
w^/L  it  varies  linearly  with  Wq/L.  This  has  very  interesting  and 
simplifying  consequences  in  the  variational  equation  as  shown  in 
Section  IVB-  The  results  obtained  by  Wes tine  and  Baker^  for  the 
rigid  - perfectly  plastic  material  are  also  shown  in  Figure  11. 


where  V denotes  an  integral  over  volume  and  6 is  the 
. . m 

max imum  s t r a x n 
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3.  Calculation  of  permanent  sets  under  impulsive  loading 


The  total  deflection (not  the  permanent  set)  of  the  beam  under  an 
impulsive  load  can  be  computed  by  the  relation  given  in  Section  IV, 

i .e. 
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P 
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2 m 


where  H 

m 

e 


m 
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P (t)  dt 
e 


f^(A)  dA 
w 


A 


(the  generalized  impulse) 

[821 


(the  generalized  mass) 


W 


P 


e 


U = energy  absorbed  ( a function  of  w^/L  ) 


The  program  given  in  Appendix  I is  used  to  perform  this  computation - 
The  results  for  the  simply  supported  beam  using  the  Bell  Theory  to 
compute  the  energy  absorption  are  shown  in  Figure  12.  The  perma- 
nent set  for  a given  impulse  is  computed  as  follows: 

a.  For  a given  value  of  impulse  in  Figure  13  pick  off  the  value 
of  Wq/L  (this  corresponds  to  the  total  deflection) 

b.  Compute  the  maximum  strain  in  the  center  of  the  beam  as  a func- 
tion of  the  center  deflection (this  curve  is  plotted  in  Fig.  14). 
At  a given  (wq/L)  center  there  is  a corresponding  value  of  (£:J^^) 
center. 

c.  Go  into  the  stress-strain  curve  of  Figure  15  with  this  value  of 
(€m)  center  and  draw  a line  parallel  to  the  elastic  line  at  this 
strain.  This  corresponds  to  the  unloading  path  from  the  maximum 
strain.  Where  the  line  hits  the  horizontal  axis  corresponds  to 
the  permanent  strain. 

d.  GO  back  to  the  Wq/L  - € curve  and  obtain  a corresponding  value 
of  Wq/L  for  this  permanent  strain.  This  will  correspond  to  the 
permanent  set  (or  permanent  deflection) . 

e.  Plot  the  impulse  parameter  originally  chosen  vs.  this  permanent 
set  parameter  as  shown  in  Fig.  12. 

5 18 

Experimental  points  obtained  from  previous  references  ' are  shown 
on  this  plot.  Indications  are  that  this  theory  gives  good  results. 
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C.  Applications  to  two  dimensional  prol)lems  - bending  and  stretching 
of  plates 

Problems  concerning  the  plastic  deformation  of  plates  involve  biax- 
ial stress  fields  and  are  therefore  an  order  of  magnitude  more  com- 
plicated to  solve  than  one  dimensional  problems.  However  the  employ- 
ment of  the  Bell  Theory  readily  enables  solutions  to  be  developed. 
The  strain  energy  for  a plate  in  a biaxial  state  of  stress  can  be 
wfitten^^ 

" • / 
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/T 


ar 


] dv 


[83] 


where  V denotes  an  integral  over  volume 


where 


T = vf  V"ox  Oy  - 3Tx^y  = 
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octahedral  shear  stress 


[84] 


= octahedral  shear  strain 

r ^ 

r = maximum  value  of  1 
m 

<J^  Oy  stress  components  , ^ ”^Y=  strain  components 

For  a plate  whose  width,  length  and  thickness  are  a,b  2H  respective- 
ly the  strain  energy  becomes 
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Now  introducing  the  Bell  Law  for  biaxial  stress,  we  have 


T = E r for  r < r 

y [86]  . 

= ^ t 5Vr-  ^ for  r .<r  ,<  r 5 = xV2  k3/2^^ 

= ^sJT  - f c ^ ^ Fc  = - ^N/X^ 

For  large  lateral  deformations  of  plates  the  strains  can  be  written 
in  terms  of  the  lateral  deflection,  w,'  as  follows: 


Now  introduce  a deflection  function  f(x,y)  such  that 

w = f(x,y) 
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where  f(x,y)  is  a nondimens ional  distribution  function.  Also  non- 
dimensionalize  the  other  parameters  as  follows: 


= x/a  , 


y = yA> 


z = 


z/H 


So  that  2the  strains  become 
= 1 ![o  8^  2 _ 2 ^ 

2 ^2  ^ 3x  a a 

€ y = 1 2 - z 

b 


The  strain  energy  -jbecomes  i 


u 


a b H E [ 


dy 


9x2 

62f 

a “2 

O y 

+1 


-1 


r = 

V w 

/ xy  = _Q  o 

a b 


r 


/ r [991 

m 


, 3f 


o 2 

z o 


— 

a b 0x  0y 


[90] 


dz 


T r,  ir  1 


[911 


Now  assume  that  rhe  edges  of  the  plate  are  simply  supported  and  let 
the  deflection  function  be  given  as ^ 


cos  [ 


^ - 1 > 
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+ cos  [* 


my  - |) 

b 

2 


J [92] 


f (x,y)  = -1(1  + cos  77-(  2 X - 1 ) ) ( 1 + cos  TT  (2  y - 1 ) ) [93] 

The  program  for  computing  the  energy  integral  is  contained  in  Appen- 
dix 2.  A plot  of  the  strain  energy  as  a function  of  the  deflection 
resulting  from  calculations  using  the  program  is  shown  in  Figure  16. 
Assuming  impulsive  loading  and  equating  the  kinetic  energy  to  the 
strain  energy  (no  corrections  are  made  for  the  shape  function  as 
done  for  the  beam  since  the  correction  is  small  for  the  simply 
supported  case)  we  can  then  obtain  a relationship  between  the  non- 
dimensional  impulse  function  and  the  nondimens ional  deflection. 

This  curve  is  shown  in  Fig.  17.  Using  an  analogous  procedure  for 
obtaining  permanent  set  as  was  used  for  beams  we  go  through  the 
following  steps: 


a.  Plot  a curve  of  nondimens ional  impulse  vs.  nondimensional  total 
deflection  such  as  shown  in  Fig.  17 

b.  Plot  a curve  of  octahedral  shear  strain,  V vs.  nondimensional 
deflection  w^/h  such  as  shown  in  Fig.  18.  For  the  given  value  of 
Wq/Yi  from  a, we  have,  a value  for  P 

c.  Go  into  the  stress-strain  curve  ( T ~ T ) in  Figure  19 

and  obtain  the  permanent  strain  by  drawing  a line  parallel  to  the 
elastic  line  at  the  strain  which  corresponds  to  the  appropriate 
value  of  w^/h. 

d.  Go  back  to  Figure  18  to  obtain  the  permanent  set  value  of  w^/h 
and  plot  on  Figure  17. 
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Both  the  nondimens ional  impulse  - deflection  values  for  total  de- 
flection and  permanent  set  are  shown  in  Figure  17  with  corresponding 
experimental  results  obtained  from  the  literature?' The  indications 
are  that  the  theory  is  quite  satisfactory  in  predicting  plastic  de- 
flections. 


D.  Bell's  theory  and  the  variational  method 


Once  the  energy  function  is  computed,  a polynomial  fit  can  be  made 
between  energy  and  deflection.  For  example,  in  the  case  of  a simply 
supported  plate  shown  in  Fig.  16  the  energy  takes  the  form 


V 

a b H E 


A (• 


+ 


B (• 


2 


[94] 


where  A = . 0009 


B = .0008 


These  coefficients  are  somewhat  less  than  those  computed  from  the 
rigid-perfectly  plastic  theory  presented  by  Westine  and  Baker.  In 
any  event  if  a second  order  polynomial  will  fit  the  energy  function 
then  all  the  theory  developed  in  section  IV  C can  be  applied  to  the 
problem  up  to  the  point  of  predicting  the  complete  isodamage  curve. 

A practical  way  to  apply  the  Bell  Theory  to  dynamic  problems  in 
plates  and  shells  is  first  to  compute  the  energy,  fit  a polynomial 
in  the  deflection  to  this  energy  function  and  then  apply  the  varia- 
tional theory  outlined  in  the  earlier  sections  of  this  report. 

E.  Possibilities  for  other  applications  of  the  Bell  theory 

The  problems  presented  in  this  section  are  only  a small  example  of 
the  application  of  Bell's  theory  to  dynamic  vulnerability  type  pro- 
blems. There  are  meiny  unsolved  vulnerability  problems  to  which  the 
new  Bell  laws  could  be  applied  to  shed  light  on  their  solution.  One 
very  important  problem  is  the  penetration  and  perforation  of  struc- 
tures by  projectiles.  These  perforation  problems  have  only  been 
investigated  either  empirically  or  analytically  using  questionable 
stress-strain  laws.  Another  problem  of  great  interest  is  the  dy- 
namic behavior  of  shells.  The  solution  of  this  problem  should  be 
merely  an  extension  of  the  analysis  presented  for  plates  in  this 
report. 
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APPENDIX  I.  COMPUTER  PROGRAM  FOR  THE  BEAM  USING  BELL'S  LAW 


A listing  of  the  beam  probram  in  "BASIC"  is  given  in  Table  1.  The 
input  parameters  are  shown  in  statements  10  - 45,  100  - 210  with 
their  corresponding  values  given  in  the  DATA  statements  480  - 600. 
The  input  parameters  in  the  program  are  given  below  with  the  sym- 
bols that  they  represent. 


C = 
El  = 
Y = 
S = 
E2  = 


= 85000  psi 
= -.285 
= .004 
= 40000  psi 
= .01 


E = elastic  modulus 


lo”^  psi 


M = number  of  integration  points  for  the  X integral  = 16 


H = number  of  integration  points  for  the  K integral  = 16 


H = h/2L  = one  half  the  thickness  of  the  beam  divided  by 
its  length  = .0069 

Wl,  W2,  W3  = minimum,  maximum,  increment  in  W where 


W = Wq/L  in  which  Wq  is  the  maximum  deflection  and 
L is  the  beam  length 

Gl,  G2,  G3  = minimum,  maximvun,  increment  in  G,  where  G is  the 
strain  value  used  to  compute  the  stress-strain 
curve  in  statements  55-97. 


X(I),G(I)  are  the  Gauss  locations  and  weights  respectively 
for  the  x integration  (locations  are  in  state- 
ments 530  - 540,  weights  in  statements  550-560) 

Y(J),  H(J)  are  Gauss  locations  and  weights  respectively  for 
the  K integration  (locations  are  in  statements 
570-580,  weights  in  590-600) 

U is ■ the  nondimensional  energy  function  v/BH'^  CJy  as  shown  in 
the  graph  of  Fig.  11 
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Table  1.  Listing  of  Beam  Program 
5 DIM  X(20)/G(205>1f(20)>H(20) 


10 

READ  C*EUY/S*E2»E 

20 

READ  M«N 

30 

READ  H 

40 

READ  U1«U2»W3 

45 

READ  C1/C2#C3 

50 

LET  Q*  0 

55 

FOR  8®G1  TO  G2  STEP  G3 

60 

IF  G>Y  THEN  72 

65 

LEI  S5®E*G 

70 

GOTO  95 

72 

IF  G>E2  THEN  90 

75 

LET  S5=S+C*SQR( G-Y) 

30 

&.TTO  9 5 

90 

LET  S5*C»SQR(G-E1) 

95 

GOTO  9? 

97 

NEXT  C 

100 

FOR  I®  0 TO  K-  1 

110 

HEAD  XCI) 

120 

NEXT  I 

130 

FOR  I®  0 TO  K-  1 

140 

READ  G(I) 

150 

NEXT  I 

160 

FOR  0®  0 TO  N-  1 

170 

READ  YCJ) 

130 

NEXT  J 

190 

FOR  o»  0 TO  N-  1 

200 

READ  H(J) 

210 

NEXT  J 

215 

FOR  V»U1  TO  W2  STEP 

W3 

220 

FOR  1=  0 TO  K-  1 

230 

FOR  0®  0 TO  N-  1 

240 

LET  F3«C  192/  5>*CX( 

I) 

250 

LET  K®H*W*F3*YC0) 

260 

IF  K>Y  THEN  290 

270 

LET  B®C  4/  3)*E*K/£ 

230 

GOTO  410 

290 

LET  A=(£/  3)*CY/K)» 

2 

300 

IF  K>E2  THEN  350 

310  LET  A1=<S/  2)*(  l-CY/K)t  g) 

320  LET  A2=<  2*C*C  2*Y+  3*K)*(K-Y)t  l.5)/( 

330  LET  B=  A>fCA  + AI+A2)/S 
3A0  GOTO  410 

350  LET  E3=(S/  25*«(E2/K)i  2-<Y/K)»  2) 

360  LET  D3*(  2*C*C  2*Y+  3*E2 ) ♦ C E2- Y ) t 1.5)/ 
370  LET  Fl=  2*C*C  2*El+  3*K)*Ol-Ent  1.5 
380  LET  F2=  2*C*(  2*El->  3*E2)  *C  E2-E1  ) » 1.5 
390  LET  F«(Fl-F2)/(  15*Kt  8) 

400  LET  B®  4*(A+B3*D3+F)/S 

410  LET  G=Q*G(I)*H( J)»F3*B 

420  NEXT  J 

430  NEXT  I 

435  LET  U«W*Q 

440  PRINT  U«U 

460  LET  e-  0 

470  NEXT  W 


15*Kt  8) 


( 15*Xf  2) 
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Table  1 - Cont'd 


430  DATA  85000# -.335 

431  DATA  4.00000E-03#  40000#  l.OOOOOE-08#  l.OOOOOE+07 

490  DATA  16#  16 

500  DATA  6.90000E-03 

5S0  DATA  5.00000E-08#  .5#  5.00000E-02 

521  DATA  l.OOOOOE-03#  l.OOOOOE-02#  l.OOOOOE-03 

530  DATA  5.29954E-03#  2.77124E-02#  6.71344E-02#  .122297 

532  DATA  .191061#  .270991#  .359193#  .452493 

535  DATA  .547506#  *640301#  .729003#  .803933 

540  DATA  .377708#  .932315#  .972237#  .9947 

550  DATA  1.35762E-08#  3.11267E-02#  4.75792E-02#  6.23144E-08 

552  DATA  7.47930E-02#  3.45732E-02#  9.13017E-02#  9.47243E-02 

555  DATA  9.47853E-08#  9.13017E-02#  3.45732E-02#  7.47930E-02 

560  DATA  6.23144E-02#  4. 75792E-02#  3.11267E-02#  1.35762E-02 

570  DATA  5.29954E-03#  2.77124E-02#  6. 71344E-02#  .122297 

572  DATA  .191061#  .270991#  .359193#  .452493 

575  DATA  .547506#  .640301#  .729003#  .308933 

530  DATA  .877702#  .932315#  .978237#  .9947 

590  DATA  1.35762E-02#  3.11267E-02#  4.75792E-02#  6.23144E-02 

592  DATA  7.47930E*02#  3.45732E-02#  9.13017E-02#  9.47253E-02 

595  DATA  9.47253E-02#  9.13017E-08#  3.45732E-08#  7.479S0E-02 

600  DATA  6.23144E-02#  4. 75792E-02#  3.11267E-08#  1.35762E-08 

700  END 
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APPENDIX  II.  COMPUTER  PROGRAM  FOR  THE  PLATE  USING  BELL'S  LAW 

A listing  of  the  plate  program  in  "BASIC"  is  given  in  Table  2.  The 
input  parameters  are  shown  in  statements  20  - 220  with  their  corre- 
sponding values  given  in  the  DATA  statements  720  - 870.  The  input 
parameters  in  the  program  are  given  below  with  the  symbols  they 


represent : 

c=  xlfR^A 

= 60000  psi 

E = elastic  modulus 

= lo”^  psi 

T1  = Ty 

= 28000  psi 

" Gl  = Ty 

= .0035 

G2  = Fc 

= .013 

G3  = Fc  - 

= -.33 

S = CTy  = 40000  psi 

M = number  of  integration  points  for  x integral  =4 

N = niamber  of  integration  points  for  y integral  = 4 

P = number  of  integration  points  for  z integral  =7 

Q = number  of  integration  points  for  F integral  = 5 

H = h/2a  = one  half  the  thickness  of  the  plate  divided  by  the 
width 

A9  = a/b  = width/length 

Wl,W2,W3  = minimiim,  maximum,  increment  in  W where  W = Wq/h 
(Wq  = maximum  deflection) 

X(I),A(I)  = locations,  weights  for  x Gauss  integration 

Y(J),B(J)  = locations,  weights  for  y Gauss  integration 

Z(K),C(K)  = locations,  weights  for  z Bauss  integration 

G(L),D(L)  = locations,  weights  for  F integration 

The  nondimensional  impulse  function  is  given  by  16  and  the  nondimen- 
sional  deflection  parameter  Wq/Ii  by  W9.  The  nondimensional  energy 
function  V/abHE  is  given  in  the  program  by  the  symbol  V. 
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Table  2.  Listing  of  Plate  Program 

10  DIM  X(  10)»y(  10)#Z(  10),G(  10)#A<  10)>BC  10)#C(  10)»DU0) 

20  READ  C.E,TUGl*G2#G3>S 

25  READ 

30  READ  H«  A9 

40  READ  kl»\«2*V3 

50  FOR  la  0 10  M-  1 

55  READ  X(I) 

60  NEXT  I 

65  FOR  I*  0 TO  K-  1 
70  READ  A<I> 

75  NEXT  I 

30  FOR  Ja  0 TO  N-  1 
35  READ  Y(0) 

90  NEXT  J 

95  FOR  0=  0 TO  N-  I 

100  READ  E(J) 

105  NEXT  0 

110  FOR  K»  0 10  R*  1 

120  READ  ZCK5 

130  NEXT  K 

140  FOR  Ka  0 TO  R-  1 

150  READ  CCK> 

160  NEXT  K 

170  FOB  L=  0 10  G-  1 

130  READ  GCL) 

190  NEXT  L 

200  FOR  L«  0 TO  G-  1 

210  READ  DCL) 

220  NEXT  L 
295  LET  V=  O 

300  FOR  WaVl  TO  W2  STEP  Vk3 

310  FOR  1=  0 TO  K-  1 

320  FOR  Ja  0 TO  N-  1 

350  LET  Xl=  3.14159*C  2»XCI)-  1) 

360  LET  Yla  3.  14159*(  24>YCJ)-  1) 

370  LET  Fl=  .25*C-6.23319*SIN(X1  J)*C  1+COSCYD) 

330  LET  F2=  .25*(-39.4734*C0S(Xl ) >♦(  1+COSCYlJ) 

390  LET  F3=  .25*(  1+COSC X 1) ) *C -6. 233 19+5INC Y 1 / ) 

400  LET  F4=  .25*C  1 ♦COSt X 1 ) ) ♦( -39. 4734+COSC Y 1 ) ) 

410  LET  F5a  .25*C-6.23319*SINCXl)>*{-6.2o319*SIN(Yl)) 

415  FOR  K=  0 TO  P-  1 

420  LEI  Ela  .5*CVit  2>»CFlt  25-ZC'AT  ♦W+F2+H 

430  LET  E2=  .5*(Gt  2)*(A9t  2)»CF3»  2)-ZCK) +H*W*A9*F4 

440  LET  E3»W*V,*A9*F1*F3-  2*ZCK> ♦H*W*A9*F5 

450  LEI  G9*  1.41399*SGR<<Elt  2)+CE2f  2)+CEl*E2)+  .25*CE3«  25) 
460  FOR  La  0 TO  G-  1 
470  LEI  G=C9*G(L) 

430  IF  G>G1  THEiN  505 

490  LET  TaG 

500  GOTO  540 

505  IF  G>G2  THEN  530 

510  LET  TaCTl/E)+CC/E)*SCRCG-Gl) 

520  GOTO  540 

530  LET  Ta(C/E)*SGR(G-G3) 

540  LET  V=U+T*G9*ACI)*B( J)*CCK)*D<L) 
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Table  2 - Cent ' d 


550  NEXT  L 
560  NEXT  X 
570  NEXT  J 
530  NEXT  I 

590  let  16-  .e5«(  l/(A9*  t2)*(£/S>>*V 

600  LET  V9-  .5*V/H 

700  PRINT  I6*W9 

710  LET  V-  0 

715  NEXT  W 

720  DATA  3S000>  1 • 00000£-»0 7«  35000«  5>00000£-03«  1*  30000£-02«  ••  33«  AOOOO 

730  DATA  A»  4#  7#  5 

740  DATA  4.00000£*02>  *59 

750  DATA  8*00000£-02*  »8*  8«00000£*02 

760  DATA  6*94313E>02«  •330009>  .66999#  .930563 

770  DATA  .173927#  .326072#  .326072#  .173927 

730  DATA  6.94318E-02#  .330009#  .66999#  .930563 

790  DATA  . 173927#  .326072#  *326072#  . 173927 

300  DATA  -.949103#-. 741531#-. 405345#  0 

310  DATA  .405345#  . .741531#  *949108 

820  DATA  .129434#  .279705#  *33133#  .417959 

830  DATA  .33133#  .279705#  .129434 

340  DATA  4.69100E-02#  .230765#  .5 

350  DATA  .769234#  .95309 

360  DATA  .113463#  .239314#  .234444 

370  DATA  .239314#  .113463 

900  END 
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SAED,  M.  Mitchell 
SAED,  D.  Tauras 
Warminster,  PA  18974 

2 Commander 

US  Naval  Air  Propulsion  Test 
Center 

ATTN:  R§T  Div,  P.  Piscopo 

ATE- 4,  Mr.  Scott 
Trenton,  NJ  08628 

1 Commander 

US  Naval  Ship  Engineering 
Center 

Prince  Georges  Plaza 
Center  Building 
Hyattsville,  MD  21402 

1 Commander 

David  W.  Taylor  Naval  Ship 

Research  5 Development  Center 
Annapolis,  MD  21402 


No.  of 

Copies  Organization 

1 Commander 

US  Naval  Surface  Weapons  Center 

ATTN:  Code  040 

Silver  Spring,  MD  20910 

4 Commander 

US  Naval  Surface  Weapons  Center 
ATTN;  Code  DEC 

Code  GAV,  T.  McCants 
W.  D.  Ralph 
J.  R.  Rech 
Dahlgren,  VA  22448 

5 Commander 

US  Naval  Weapons  Center 
ATTN:  Code  408,  H.  Drake 

Code  4070,  M.  Keith 
Code  4081,  C.  Sandberg 
Code  6031,  M.  Backman 
Code  4083,  G.  Monscko 
China  Lake,  CA  93555 

1 Commander 

US  Naval  Ammunition  Depot 
ATTN:  RD-3 

Crane,  IN  47522 

1 Director 

US  Naval  Research  Laboratory 
ATTN:  Code  5554,  Dr.  Esterowitz 

Washington,  DC  20390 

3 Commandant 

US  Marine  Corps 
ATTN : Code  AAP 

AAW-1,  LTC  Peterson 
AAW-5,  LTC  Paige,  Jr. 
Washington,  DC  20380 

1 Director 

Development  Center,  MCDEC 
Quant i CO,  VA  22134 

1 HQ  USAF  (AFSCAGF) 

Washington,  DC  20330 
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3 AFSC  (LTC  Stephenson) 

Andrews  AFB 
Washington,  DC  20331 

1 AFSC  (SDOA) 

Andrews  AFB 
Washington,  DC  20331 

3 HQ  USAFSC  (XRLW,  LTC  Seufert; 
SDW,  R,  Hartmeyer;  DLCA, 

P.  Sandler) 

Andrews  AFB 
Washington,  DC  20334 

1 USAFTAWC  COA) 

Eglin  AFB,  FL  32542 

3 AFATL  (DLY) 

Eglin  AFB,  FL  32542 

1 AFATL  (DLRW,  W.  D.  Thomas) 

Eglin  AFB,  FL  32542 

1 OSU 

ATTN:  Mr.  R.  Armstrong 

Eglin  AFB,  FL  32542 

2 AFWL  (MAJ  H.  Rede; 

MAJ  J.  Godsey) 

Kirtland  AFB,  NM  87117 

2 AFSWC  (TEEM,  LTC  D.  Meyers; 

TESPO,  MAJ  R.  Gerhandt) 
Kirtland  AFB,  NM  87117 

1 TAG  (DIO) 

Langley  AFB,  VA  23365 

4 ASD  (ASBR;  ASBO,  Mr.  Erkeneff; 

ASZT;  XOO,  Mr.  Burneka) 
Wright-Patterson  AFB,  OH  45453 

3 ASD  (ENYS,  J.  Wallick; 

J.  Howard;  L.  DiRito) 
Wright-Patterson  AFB,  OH  45433 


No.  of 

Copies  Organization 

1 ASD  (XRHL,  G.  B.  Bennett) 
Wright-Patterson  AFB,  OH  45433 

4 AFLC  (MCNEA,  W.  Dougherty; 

PPSC,  R.  Mergler;  MMHM, 

E.  Swanson;  MEA,  COL  T.  Egan) 
Wright-Patterson  AFB,  OH  45433 

4 AFAPL  (SFL,  Jones;  SFH,  Botteri; 
SSF,  Churchill;  AFH, 

R.  Clodfelter) 

Wright-Patterson  AFB,  OH  45433 

5 AFFDL  (FBS;  FER;  PTS,  D.  Voyls; 
R.  Lauzze;  CPT  Beardon) 

Wright-Patterson  AFB,  OH  45433 

2 AFML  (MAE;  MAS,  COL  Hughes) 
Wright-Patterson  AFB,  OH  45433 

2 Director 

National  Aeronautics  and 
Space  Administration 
Langley  Research  Center 
ATTN:  MS  246E,  L.  Bement 

MS  249,  J.  Ward 
Hampton,  VA  23665 

1 Director 

National  Aeronautics  and 
Space  Administration 
Lewis  Research  Center 
Lewis  Directorate 
Cleveland,  OH  44135 

1 Bell  Helicopter  Company 
ATTN:  Mr.  Nile  Fischer 

P.O.  Box  482 
Fort  Worth,  TX  76101 

1 The  Boeing  Company 
Vertol  Division 
ATTN;  Dave  Harding 
P.O.  Box  16858 
Philadelphia,  PA  19142 
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1 Fairchild  Industries 

Fairchild  Republic  Division 
ATTN:  D.  Watson 

Farmingdale,  LI,  NY  11735 

1 Falcon  Research  and 

Development  Company 
ATTN:  Lavelle  Manhood 

1225  South  Huron 
Denver,  CO  80223 

1 Falcon  Research  and 

Development  Company 
ATTN:  Arthur  Stein 

One  American  Drive 
Buffalo,  NY  14255 

1 General  Electric  Company 

Military  Engine  Products  Div 
ATTN:  E.  L.  Richardson 

Building  2406A 
1000  Western  Avenue 
West  Lynn,  MA  01910 

1 Summa  Corporation 

Hughes  Helicopter  Division 
ATTN:  J.  H.  Tinley 

Centinela  and  Teale  Streets 
Culver  City,  CA  90230 


No.  of 

Copies  Organization 

1 LTV  Aerospace  Corporation 
Vought  Systems  Division 
ATTN:  D.  M.  Reedy  - Unit  2 

54244 

P.O.  Box  5907 
Dallas,  TX  75222 

1 Rockwell  International 

Los  Angeles  Aircraft  Division 
ATTN:  R.  L.  Moonan 

Mail  Code  AB75 
International  Airport 
Los  Angeles,  CA  90009 

1 United  Aircraft  Corporation 
Sikorsky  Aircraft  Division 
ATTN:  J.  B.  Foulk 

Stratford,  CT  06602 


Aberdeen  Proving  Ground 

Marine  Corps  Ln  Ofc 
Cdr,  USATECOM 

ATTN:  Mr.  Moshang 

Dir,  USAMSAA 

ATTN ; J . Lindenmuth 
D.  O'Neill 
J . McCarthy 
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